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Abstract 

We study solutions of the 2D Ginzburg-Landau equation 



vn ■ -Au + ^n(|np-l) =0 

o 



subject to " semi-stifF' boundary conditions: the Dirichlet condition for the 
modulus, \u\ = 1, and the homogeneous Neumann condition for the phase. 



. The principal result of this work shows there are stable solutions of this pro- 

blem with zeros (vortices), which are located near the boundary and have 



bounded energy in the limit of small e. For the Dirichlet bondary condi- 



^ I tion ("stiff' problem), the existence of stable solutions with vortices, whose 

energy blows up as e — > 0, is well known. By contrast, stable solutions with 
vortices are not established in the case of the homogeneous Neumann ("soft") 
boundary condition. 

In this work, we develop a variational method which allows one to construct 
local minimizers of the corresponding Ginzburg-Landau energy functional. We 
introduce an approximate bulk degree as the key ingredient of this method, 
and, unlike the standard degree over the curve, it is preserved in the weak 
iJi-limit. 
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1 Introduction and Main Results 



In this work, we study solutions of the Ginzburg-Landau (GL) equation 

-Au + \u(\u\^-l)^0 in A, (1.1) 

where £ is a positive parameter (the inverse of the GL parameter /t = 1/e), m is a 
complex-valued (R^-valued) map, and A is a smooth, bounded, multiply connected 
domain in M^. For simplicity, hereafter we assume A is an annular type (doubly 
connected) domain of the form A = fl\uj, where fl and uj are simply connected 
smooth domains and U C ^2 C M^. 

Equation (1.1) is the Euler-Lagrange PDE corresponding to the energy functional 

Eeiu) = I ^ I V^l^o; +±J^{\u\'- Ifdx. (1.2) 

Equations of this type arise, e.g., in models of superconductivity and superfluidity. 
Additionally, (1.1) is viewed as a complex- valued version of the AUen-Cahn model 
for phase transitions [32] . 

Solutions of (1.1) subject to the Dirichlet boundary condition, u — g on dA 
with fixed (S^-valued boundary data g, have been extensively studied in the past 
decade. Special attention has been paid to solutions with isolated zeros (vortices). 
In contrast with the Dirichlet problem, in the case of the homogeneous Neumann 
boundary condition, solutions are typically vortexless; in particular, stable solutions 
with vortices have not been established. 

This work is devoted to solutions of (1.1) subject to the "semi-stifF' boundary 
conditions 

du 

\u\ = 1 and w X — = on dA. (1.3) 

These boundary conditions are intermediate between Dirichlet and Neumann in the 
following sense: any solution u G H^{A;W'^) of (1.1, 1.3) is sufficiently regular [11], 
so it can be written as u = \u\&^ (locally) near the boundary. Then (1.3) means 
the Dirichlet boundary condition is prescribed for the modulus, \u\ = 1 on dA, and 
the Neumann condition is prescribed for the phase, |^ = on dA. 

Problem (1.1, 1.3) is equivalent to finding critical points of the energy functional 
(1.2) in the space 

J = {ue H\A;R^y, \u\ = 1 a.e. on OA}. (1.4) 
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Our main objective is to study the existence of stable solutions of (1.1, 1.3) with 
vortices. Since the problem is time independent, stable solutions are defined as 
(local) minimizers of (1.2) in J'. In other words, we are interested in whether the 
model (1.1, 1.3) stabilizes vortices similarly to Dirichlet problem or does not stabilize 
vortices analogously to Neumann problem. The boundary conditions (1.3) are not 
well studied, and this work, along with studies [7,11-13,21] reveals their distinct 
features, described later in the introduction. 

Let us briefly review the existing results for the Dirichlet and Neumann boundary 
value problems for equation (1.1). The first results on the existence of stable solu- 
tions with vortices for Dirichlet problem were obtained in [19,20]. Stable solutions 
of (1.1) with vortices were obtained and studied in [14] for star-shaped domains and 
prescribed -S'^-valued boundary data with nonzero topological degree. In [14], the 
limiting locations (as £ — > 0) of vortices of globally minimizing and other solutions 
(if they exist) are described by means of a renormalized energy. Subsequently, these 
results were generalized for multiply connected domains in [34]. The existence of 
locally minimizing and minmax solutions was established first in [28] and [29] , then 
in more detail and generality in [27] (see also [4,17]). We refer the reader to [9], 
and references therein, for the Dirichlet problem's various results. As previously 
mentioned, only vortexless stable solutions of (1.1) with the homogeneous Neumann 
boundary condition in 2D are known. Moreover, all locally minimizing solutions 
arc constant maps if A is convex [22], or simply connected and e is small [33]. The 
existence of nonconstant (but vortexless) locally minimizing solutions is established 
in [23] and [3]. In the recent work [17], a general result for the existence of (nonmi- 
nimizing) solutions with vortices was obtained. Similarly to the Dirichlet problem, 
these solutions with vortices have energy that blows up as e ^ 0. 

Equation (1.1) (functional (1.2)) is usually referred to as a simplified GL model 
(without magnetic field). There is a large body of mathematical literature on the 
general GL model with a magnetic field (e.g., [2,5,24,30,31]). Since (1.1) is obtai- 
ned from the general GL energy by setting the magnetic field to zero, it describes 
persistent currents in a 2D cross-section of a cylindrical superconductor (or in a 
2D film). It was observed in [14], the degree of the boundary data on connected 
components of dA creates the same type of " quantized vortices" as a magnetic field 
in type II superconductors or as angular rotation in superfluids. Despite a relatively 
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simple form of equation (1.1), it leads to a deep analysis of properties of its solutions 
similar to other fundamental PDE's in mathematical physics. 

The boundary conditions (1.3) model, e.g., the surface of a superconductor coa- 
ted with a high temperature superconducting thin film [6]. Generating a mathema- 
tical model of persistent currents in such a superconductor, then amounts to finding 
critical points of functional (1.2) in the space J ^ when u—\u\e^'^ow. the boundary 
and \u\ = 1, while the phase ip is "free". 

Boundary conditions (1.3) appeared in recent studies [7,12,21] of the minimi- 
zation problem for the GL functional (1.2) among maps from J' with prescribed 
degrees on the connected components of the boundary. The minimization of the 
energy (1.2) in J produces only constant solutions of (1.1, 1.3), similar to the case 
of the Neumann problem, which corresponds to finding critical points of (1.2) in the 
entire space if^(A;]R^). An obvious way of producing critical points with vortices 
is to impose two different degrees q p on dfl and du. That is, to consider the 
minimization of Es{u) in the set C JT", where 



where the integral is understood via H^^"^ — H~^/'^ duality, and ^ is the tangential 
derivative with respect to the counterclockwise orientation of 7. (Throughout the 
paper we assume the same orientation of duj and dVt.) Note that Jpq are connected 
components of J (see [8]). 

Simple topological considerations imply that critical points from Jpq must have 
at least |p — q'| (with multiplicity) vortices. We emphasize that the existence of 
such critical points is far from obvious. For example (see Section 2), there are no 
global minimizers of E^{u) in jToi and the weak limits of minimizing sequences do 
not belong to jToi. This simple example illustrates an important property of the sets 
Jpq, which is crucial for our consideration: these sets are not weakly -closed, since 
the degree at the boundary may change in the limit. On the other hand, the results 



Jpq {u e J] deg(ii, duj) = p, deg('u, dVt) = q.} 



(1.5) 



Recall that the degree (winding number) of a map u G if^/^(7, S^) on 7 (where 7 is 
either duj or dVL) is an integer given by the classical formula (cf., e.g., [8]) 




(1.6) 
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of this work show that when p = q and there is no topological reason for vortices to 
appear, local minimizers typically do have vortices. 

As mentioned above, the vortex structure of solutions of (1.1) with Dirichlet 
and Neumann boundary conditions is well studied. In contrast, only vortexless 
solutions of the semi-stiff problem (1.1), (1.3) were found [21], [12]. In [12], it 
was shown that minimizing sequences for the corresponding minimization problem 
develop a novel type of so-called "near-boundary" vortices, which approach the 
boundary and have finite GL energy in the limit of small e (due to the ghost vortices, 
see Appendix A). However, such minimizing sequences do not converge to actual 
minimizers [13]. These studies lead to the natural question of whether there exist 
true solutions of (1.1, 1.3) with near-boundary vortices. Unlike the minimizing 
sequences such solutions may model observable states of a physical system (e.g., 
persistent currents with vortices and superfluids between rotating cylinders [18]). 
The following theorem, which is the main result of this work, provides the answer 
to this question. 

Theorem 1 (Existence of solutions with vortices of problem (1.1, 1.3)). For any 

integer M > 0, there exist at least M distinct stable solutions of (1.1, 1.3) with 
(nearhoundary) vortices when e < Si (si — £i(M) > Q). The vortices of these 
solutions are at distance o{e) from the boundary and have hounded GL energy in the 
limit £ — > 0. The solutions are stable in the sense that they are (local) minimizers 



To construct local minimizers of (1.2) in J , we represent J as the union of 
subsets Jp'^ (defined in (1.8) below), J — Up^q^dez>Jpq\ and study the existence of 
global minimizers in J'pq^ ■ Furthermore, we show that each minimizer lies in J'pg^ 
with its open neighborhood. Therefore, the minimizers in J'pq^ are distinct local 
minimizers in J^. 

Thus, the construction of solutions of (1.1, 1.3) is based on the study of the 
following constrained minimization problem: 



of (1.2) in J. 



me{p, q, d) iTd{E^{u)] u G J^f>}, 



(1.7) 



where 



J^f = {m e Jpq] abdeg(M) e [rf - 1/2, d + 1/2]}, 



(1.8) 
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p,g and d are given integers, and abdeg(-u) is the approximate hulk degree, introduced 
as follows. Consider the boundary value problem 

AT/ = in A 
< V =1 ondQ (1.9) 
V = Q on duj. 
Introduce abdeg( • ) : H^{A;R'^) ^ M by the formula 

abdeg(ti) = [ u x {d^^V d^^u - d^^V d^^u) dx, (1.10) 

^TT J A 

where V solves (1.9). In the particular case where A is a circular annulus, Ar^r^ — 
{x;Ri < \x\ < R2}, abdeg(M) is expressed by 

abdeg(ti) = , / , , [""'(J- f ux ^ds) ^. (1.11) 
' log(i?2/i?i) V27r4|=^ dr J ^ ^ ' 

For S'^-valued maps, abdeg(w) becomes integer valued and representation (1.11) 
clarifies its interpretation as an average value of the standard degree. The definition 
(1.10) is motivated by the following intuitive consideration: represent the standard 
degree over the boundary dQ via a " bulk" integral over the area of A for S'^-valued 
maps and notice that if E^{u) < A for some finite A and sufficiently small e, then u 
is " almost" 5" valued. 

It was observed in [3] that for 5'^-valued maps in an annulus A, one can define the 
topological degree deg{u, A) that classifies maps u e H^{A; S^) according to their 1- 
homotopy type [35] (1-homotopy type is completely determined by the degree of the 
restriction to a nontrivial contour) . This definition was relaxed in [3] for maps that 
are not necessarily 5" valued by considering u/\u\ in a subdomain A^ C A, which is 
obtained by removing neighborhoods of the boundary dA and zeros (vortices) of u. 
Definition (1.11) does not require the removal of vortices from A, and abdeg{u) is 
obtained by a simple formula (unlike deg(M, A) in [3] , where the domain of integration 
depends on u). Note that in general abdeg(M) is not an integer. The most important 
fact for our consideration is that abdeg(M) is continuous with respect to weak H^- 
convergence, unhke the standard degree in (1.6) (this issue for deg{u,A) was not 
addressed in [3]). 

The minimization in problem (1.7) is taken over J^^q\ which is not an open set, 
and therefore minimizers of (1.7) (if they exist) are not necessarily local minimizers 
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of (1.2) in J . Indeed, while J^q is an open subset of J (hereafter we assume the 
topology and convergence in J to be the strong unless otherwise is specified), 
the constraint abdeg(-u) G [rf — 1/2, + 1/2] defines a closed set with respect to both 
strong and weak i7^-convergences. However, if we further consider a subset of maps 
with bounded energy and choose e small enough, then the constraint abdeg(M) € 
[d — 1/2, 0? + 1/2] becomes open thanks to the following proposition: 

Proposition 2. Fix A > 0. There exists Eq = eo{A) > such that if < e < So, 
then for any integer d and any u e H^{A;M.'^) satisfying Es{u) < A the closed 
constraint abdcg(u) E [d — 1/2, d+ 1/2] is equivalent to an open one that is, 

d-l/2<abdeg{u) <d+l/2 <s=^ d - 1/2 < abdeg{u) < d + 1/2. (1.12) 

The following theorem is the main tool in proving existence of local minimizers. 

Theorem 3 (Existence of minimizers of the constrained problem). For any integers 
p, q and d > (d < 0) with d > max{p, g} (d < mm{p,q}) there exists £i = 
ei{p,q,d) > such that the infimura in (1.7) is always attained, when e < e^. 
Moreover 

me{p,q,d) < Io{d, A) + 7r{\d - p\ + \d-q\), (1.13) 

where 

Io{d,A)^mmi^^ JjWu\^dx, u e H\A, S') D (1.14) 

The value Io{d,A) is expressed by Io{d,A) = 2{7idy / ca.p{A) via the H^-capacity 
cap(yl) of the domain A. 

The key difficulty is to estabhsh the attainability of the infimum in (1.7), which 
is highly nontrivial since the degree on and duj is not preserved in the weak H^- 
limit [12,13]. We show solutions of (1.1, 1.3), which are minimizers of (1.7) (local 
minimizers of (1.2) in J) with p ^ d and any q {pi q ^ d and any p) must have 
vortices. For fixed £, these vortices are located at a positive distance from dA and 
approach dA as £ ^ 0. 

Without loss of generality, throughout this work we always assume that d > 
(otherwise one can reverse the orientation of M?). 

Theorem 1 follows from Theorem 3 and Proposition 2. The asymptotic behavior 
of the local minimizers is established in 
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Theorem 4 (Asymptotic behavior of minimizers and their energies). Assume that 
the integers p, q and d satisfy the assumptions of Theorem 3. Then as e ^ Q 
minimizers of (1.7) converge weakly in H^{A), up to a subsequence, to a harmonic 
map u which minimizes (1.14). Additionally, 

Es{ue) — Io{d, A) + 7r{\d — p\ + \d — q\) + o{l), as e ^ 0, and (1-15) 

Eeiue)^}: \Vu^\'^dx + o{l), as E ^ 0. (1.16) 
2 J A 

In particular, it follows from (1.15, 1.16) that there is no strong convergence of 
minimizers of (1.7) in H^{A) as £ ^ unless p = q = d. 

Next, we summarize the distinct features of the GL boundary value problem 
with semi-stiff boundary conditions. The first interesting feature is the existence of 
solutions with a new type of vortices called near-boundary vortices. Unlike the inner 
vortices, whose energy blows up at the rate of | log£|, the energy of near-boundary 
vortices is bounded as £ ^ and they are located at a distance o{e) from the 
boundary. 

Secondly, the semi-stiff boundary conditions result in a lack of compactness. 
Namely, as of now, the only way to find nonconstant minimizers is by searching for 
minimizers in subsets C J. These subsets, however, are not weakly iJ^-closed 
and therefore a weak if^-limit of a minimizing sequence (u^''^) e J^^ may not lie in 
but rather in J^f^, with p' ^ p or (and) q' ^p. Theorem 3 shows that if d > 
and d > max{p, q}, then (for small e) any weak if^-limit of a minimizing sequence 
{u^'^^) C Jp^ always belongs to Jp^ , despite the lack of weak if ^-closeness of Jpf. 
In contrast, if d > and d < max{p, q} we have 

Conjecture 5. Let d > 0, d < max{p, q} {or d < 0, d > min{p, q}) and let m be a 
weak limit of a minimizing sequence for problem (1.7) (such a minimizing sequence 
exists and bound (1.13) holds for any integer p, q, d, see Appendix B). Then u ^ Jp^ 
when £ is sufficiently small. 

In the simplest case, when d = and either p= \ and q — ^orp — ^ and g = 1, 
this conjecture is demonstrated by an argument quite similar to the nonexistence 
proof in [11] for simply connected domains (see Sec. 2 below). A more interesting 
example, which supports the above conjecture, follows from the previously studied 
(global) minimization problem = mf{E^{u), u G Jn}. It was shown in [11,13] 
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that if cap(A) > tt (subcritical/critical cases), then is always attained, whereas 
if cap(A) < TT (supercritical case), then is never attained for small e. One can 
see elements of minimizing sequences lie in J'H^ in subcritical/critical cases and in 
J'l'l^ in the supercritical case. Moreover, the nonexistence of minimizers in problem 
(1.7) for d = (},p = q = l and small e holds for any doubly connected domain (with 
any capacity). (For cap(A) < tt the proof is presented in [13], this proof can be 
easily generalized for ca.p{A) > vr.) 

We conclude the introduction by outlining the scheme of the proof of Theorem 3, 
which employs a comparison argument. Fix an integer d > 0. First, we establish the 
existence of minimizers in problem (1.7) ior p — q — dhy using the so-called Price 
Lemma [12] (see Lemma 9 below), the uniform lower energy bound from Lemma 16 
and the upper bound from Lemma 14, which is obtained by considering 5'^-valued 
testing maps. We show these minimizers (which belong to J^jf) are vortexless. 
Next, we argue by induction on the parameter ae(p, g) = \d — p\ + \d — q\. This 
parameter is naturally associated with the number of vortices- for example, for the 
above minimizers in J^^f, we have ee(d, d) — 0. Given an integer X > 0, we assume 
the existence of minimizers in problem (1.7) for p, q such that, ae(p, g) < K and 
P ^ d, q < d (the induction hypothesis) and prove the existence of minimizers for 
p, q such that ae(p, q) = K + 1 and p < d, q < d. The first step in the induction 
procedure (when K = 0) is shown in Section 5. The key technical point there is to 
construct a testing map v G J^^i^d-i) ^^^^ 

E,{v)<E,{uo)+n, (1.17) 

where uq is a minimizer of (1.2) in J^jf- This map v is constructed by using the 
minimizer uq and Mobius conformal maps (Blashke factor [15]) on the unit disk 
with a prescribed single zero near the boundary. Then, given a minimizing sequence 
(u^'^)) C J^^^_i) of problem (1.7) ior p — d, q — d — 1, we have, by (2.7) from Lemma 
[9] and (1.17), 

E^{u)+TT{\d-deg{u,du;)\ + \d-l-deg{u,dn)\) < lim E^iu^''^) < E^{uo)+tt, (1.18) 

fe— >oo 

where u is a weak //-"^-limit of (u'^''^) (possibly a subsequence). Then we estimate the 
left hand side of (1-18) by the lower energy bound from Lemma 16 and the right 
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hand side of (1-18) by the upper bound from Lemma 14, this yields 

Io{d,A)+Tr{2\d-deg{u,duj)\ + \d-l-deg{u,dft)\ + \d-deg{u,dn)\) < Io{d,A) + ^Tr. 

(1.19) 

This imphes that deg{u, du) — d, and either deg{u, dVt) — d—1 or deg{u, dQ) — d. In 
view of (1.18), the only possible case is actually deg(M, dQ) = d—1 since, otherwise, 
u e J^f and therefore Ee{u) > Ee{uo) which contradicts (1.18). Thus u e 
and is a minimizer in J'^^}_iy The proof of existence of minimizers for p = d — 1, 
q = d is quite similar. So wc have shown that the existence of minimizers for 
3s(p, Q') = implies the existence of minimizers for ae(p, g) = 1. In the general case, 
when passing from ae(p, q) < K to a2{p, q) < K-\-l in problem (1.7), wc use the same 
idea but it is technically much more involved. It requires the asymptotic analysis as 
£ — * of minimizers Upq of (1.7) with ae(j9, q) = K, which is carried out in Section 6. 
Based on the result of this asymptotic analysis, wc construct testing maps v G J'^f^, 
ip' = p, q' = q — 1 or p' — p, q' — q — 1), such that E^{v) < E^{up^ + tt. 



2 Preliminaries 

Throughout the paper we use the following notations. 

• The vectors a = (ai, 02) are identified with complex numbers a = ai + ia2. 

• a ■ b stands for the scalar product a ■ b — aibi + 02^2 = |(o6 + ab). 

• a X b stands for the vector product a x b = 0162 — 02^1 = ^{o,b — db). 

• The orientation of simple (without self intersecting) curves in (in particular 
duj and dO) is assumed counterclockwise. If C is such a curve, r stands for the 
unit tangent vector pointing in the sense of the above mentioned orientation 
on £, u is the unit normal vector such that (z/, r) is direct. 

• If /i is a scalar function, then V-^h — {—dx^h, dx^h). 

2.1 Properties of solutions from problem (1.1, 1.3) 

As shown in [11] by a bootstrap argument, any solution u G H^{A\M?) of problem 
(1.1, 1.3) is sufficiently regular (e.g., u e C^(/l) if A has a boundary). By the 
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maximum principle we also have 

Lemma 6. The function p{x) — \u{x)\ satisfies p < 1 in A. 

Locally, away from its zeros, u can be written as u = pe^'^ with real- valued phase 
4>. We also will frequently make use of the current potential h related to the solution 
u of (1.1, 1.3) by 

W^h — (u X dx,u, u X dxou) in A 

(2.1) 

h — 1 on dQ. 

Unlike the phase (p, the function h is defined globally on A, and 

V/i = -p^V^cf) when p > 0. (2.2) 

The existence of the unique solution of system (2.1) and its elementary properties 
are established in the following 

Lemma 7. There exists the unique solution h of the system (2.1) and h = Const 
on duj, moreover 

Ah = 2dx^u X dx^u in A, (2.3) 
div(-^V/i) = when p > 0. (2.4) 

Proof. The vector field F = {u x dx^u,u x dx^u) is divergence free. Indeed, 
since w is a smooth solution of (1.1), we have divF = u x An = in A. It follows 
that for any simply connected domain W G A there is a unique (up to an additive 
constant) function $ solving V-'"$ — F inW (this is well known Poincare's lemma) 
. Such a local solution $ can be extended to a (possibly multi-valued) solution on 
A. Thanks to the fact that u satisfies (1.3) we have |^ = —F • u — on dA, i.e. $ 
takes constant values on every connected component of the boundary. This means 
$ is actually a single valued function. Then, h{x) — $(x) — $(9Q) -|- 1 is the unique 
solution of (2.1). 

The verification of (2.3) is straightforward; (2.4) follows directly from (2.2). □ 
In what follows, we also use the following result, which is valid for any solution 
of the GL equation (1.1) (not necessarily satisfying (1.3)). 
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Lemma 8 ([26]). Let u be a solutions of the GL equation (1-1) such that \u\ < 1, 
Eeiu) < A, where A is independent of e. Then 



and 

I^V^)I<— (2.6) 
dist (x, a A) 



where C, are independent of e. 

2.2 Minimization among maps of J with prescribed degrees 

Any minimizer of (1.2) over the set Jpq with prescribed integer degrees p and q is 
clearly a solution of (1.1), (1.3). However, the existence of minimizers is a nontrivial 
problem. In [7, 11-13,21] the minimization problem for the Ginzburg-Landau func- 
tional (1.2) in jTii was considered. In the case when A is a circular annulus, it was 
observed in [7] that minimizers, if they exist, brake the symmetry when the ratio of 
the outer and inner radii of the annulus exceeds certain threshold. By contrast, in 
the case when this ratio is sufficiently close to 1, the existence of a unique minimizer 
and its symmetry is shown in [21]. The techniques in both [7] and [21] rehed on 
the circular symmetry of the domain. A more general approach based on the Price 
Lemma was proposed in [11], [12]. 

Lemma 9. [12] Let {u^''^ e J'pq) be a sequence that converges to u weakly in 
H^{A,m?). Then 

liminfi / \Vu^^'>\^dx >- [ \Vu\^dx + n{\p - dcg{u,du)\ + \q - dcg{u,dQ)\), 
^ Ja Ja 

or, equivalently (by Sobolev embeddings) , 

\mi\ni E^{u^^^) >E,{u) + T:{\p-deg{u,duj)\ + \q-deg{u,da)\). (2.7) 

k 

This result is also of prime importance in this work. With the help of Lemma 
9, it was shown in [12] that the infimum of (1.2) in jTii is always attained when 
cap (A) > TT. It was also conjectured in [12] that when cap (A) < tt and e is sufficiently 
small the weak limit of any minimizing sequence is not in the class of admissible 
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maps, i.e. the global minimizer does not exist. In [13] this nonexistence conjecture 
was proved by a contradiction argument based on explicit energy bounds. 

While the existence/nonexistence of minimizers in J^q for p = g = 1 is nontrivial 
and the answer depends on cap(74) and also on e, the case p = 0, g = 1 (p = 1, 
g = 0) is simple. Arguing as in [11] we can show that \vl{Es{u)\u e jToi} is never 
attained. Really, we have 

— 7r|deg(M, 50) — deg(M, diS) \ — n 

whenever u G Joi- On the other hand, by constructing explicit sequence in the 
spirit of [7] (see also [11]), we have mi{Es{u);u G Joi} = ^r. Thus, if there exists a 
minimizer u G Joi then u G H^{A;S^) and u solves the GL equation (1.1). This is 
impossible unless w is a constant map, and then u ^ Jq\. 



IVwrdx > 



3 Properties of the approximate bulk degree 

The degree of restriction of maps from H.^^A^S^') to any smooth closed curve, in 
particular is preserved by the weak iJ^-convergence. This result follows from 
[35], or can be shown directly by using integration by parts as below in (3.2) (note 
that, for any jS'^-valued map li, deg(u, 9fi) = deg(ii, dio) — deg(ii, £), where £ is an 
arbitrary smooth simple curve in A enclosing a;). Thus we have the decomposition 

E^{A,S^) = \^{u^E^{A,S^), deg(M,aO)=d} (3.1) 

by disjoint sets, each of them being closed in weak iJ^-topology of i/^(A, 

Fix A > 0. In this section we consider maps u G i/^(A, R^) in the level set 
= {u; E^{u) < A}. We show that the approximate bulk degree abdeg(-u) classifies 

maps u G E^ similarly to the above classification (3.1) of 5'^-valued maps. The basic 

properties of abdeg(ti) we demonstrate are 

a) abdeg(ii. A) = deg(u, dQ) Hue H\A, S^), 

b) |abdeg(-u) — abdeg(f )| < |||\/||ci(yi)A^/^||-u — f 11^2(^1) if u,v e E^. 

The first property follows directly from the definition (1.10) of abdeg(M). Really, 

integrating by parts in (1.10), we get 

If du If 

abdeg(M) = — / ux -^ds / d^j^u x d^^uVdx — deg{u,dfl), (3.2) 

27r Jan or tt 
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for any u G H^{A, S^) {dx^u x dx^u = a.e. in A since \u\ = 1 a.e.). The property 
b) of abdeg(M) is proved in the following 

Lemma 10. For any u,v e {{^{A-jM."^) we have 

|abdeg(x.) - abdeg(^)| < + - ^lUv)- 

TT 

Proof. We have, integrating by parts, 
27r(abdeg(M) — abdeg(t')) = J {u — v) x {dx2udx^V — dx^udx2V)dx 



+ vx {d^^iu - v) d^,V - d^,{u - v) d^^V)^^^ 

J A 

= / (u-v) X {dx^udx^V - dxj^udx2V)dx 
J A 

+ J (u-v) X {d^^vdx^V - dx^vdx2V)dx. 



Then the statement of the lemma follows by the Cauchy-Schwartz inequality. □ 
The main consequence of properties a) and b) of the function abdeg(M) is 

Proposition 11. abdeg(w) is close to integers uniformly in u E when e is 
sufficiently small, i.e. 

c) sup dist(abdeg(M), Z) — >• as £ — >• 0. 

Before proving this fact note that Proposition 11 immediately implies Proposition 
2 stated in the Introduction. 

Proof of Proposition 11. According to (3.2) and Lemma 10 we have 

2 1 

sup dist(abdeg(M), Z) < sup inf |abdeg(M) — abdeg(u)| < — ||V^||ciA25g, (3.3) 

where is the following (nonsymmetric) distance 

5e := sup distL2(A){u,EQ) (3.4) 



between E!^ and 



E^ = |m e H\A,S^y, Eo{u) ^^J \Vu\Mx < a| . 
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Show now that 6s ^ as e ^ 0. In view of (3.3) this yields the desired result. 
Assume by contradiction that ^e^. > c > for a sequence Sk ^ 0. By virtue 
of the Sobolev embeddings the supremum in (3.4) is attained on E^, i.e. 5s = 
dist ^2 (^j^-j {us, Eq), where E^iu^) < A. We can extract a subsequence of {usj, still 
denoted (m^^), that converges to a map u weakly in H^{A,M.'^). Thanks to Sobolev 
embeddings — > u strongly in L'^{A) and we have u G H^{A, S^), since Xid^ep — 
l)^dx < 4Ae^ . Besides Eq[u) < A, by the lower weak semicontinuity of the Dirichlet 
integral. Thus u G E^ and 5^^, < \\u — Us^.\\l'^(^a) —^0. □ 
The following result illustrates the relation of abdeg(w) , which is not necessarily 
an integer with the standard notion of degree on a curve. It provides a simple 
criterion whether the constraint abdeg(w) G [d— 1/2, d-\- 1/2] in (1.7) is satisfied in 
a particular case when u is a solution of equation (1.1). 

Lemma 12. Let C = {x E A; V{x) = 1/2} denote the 1/2 level set of V , where V 
is the solution of (1-9). (C is a smooth curve enclosing oo.) Then if a solution u of 
the GL equation (1.1) satisfies \u\ < 1 in A and Es{u) < A, then (i) \u\ > 1/2 on C 
and (ii) we have 

11 

abdeg(M) e [d - 1/2, d + 1/2] deg(^, £) = d, (3.5) 

when e < e\, where Si — £i(A) > does not depend on u. 
Proof. Consider the domain 

As. = {xeA- S' < V{x) < 1 - 6'}, (3.6) 

where < 5' < 1/2. It follows from Lemma 8 that u satisfies 

\u\ > 1/2 on As', (3.7) 

when e < e[ {e'l = e[{S') > 0). This proves (i). We can now write u — pe*^ 
(p = \u\ > 1/2) on As'. Find an extension of t/j onto the whole domain A. To this 
end we pass to a conformal image of A. 

It is well known (see,e.g. [1]) that there is a conformal mapping Q of A onto 
the annulus O with R = exp(7r/cap(A)) and 1/i? as the outer and inner radii, 
correspondingly. Moreover, Q is explicitly given hj Q = exp{-^^^^({V — l/2 + i\E')), 
where ^' is a (multivalued) harmonic conjugate of V. Q maps A^' onto the annulus 
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Q{As') C O whose outer and inner radii are R' = exp(^^|py(l/2 — 6')) and l/R', 
correspondingly. 

Now consider ^(x) = ^(G~^{x)) on Q{As'). Wc can extend ip to the whole O by 
reflections iJj{x) := when > R' and ip{x) := ip{x/{R'\x\y) when 

\x\ < l/R', so that 

[ |VV'|Ma;< [ |VV'pdx+ [ |V^|Mx < 2 /" |V^|Ma;, (3.8) 

when < 5' < 1/4. Then (3.8), the conformal invariance of the Dirichlet integral 
and (3.7) imply 



L 



\Vi^\^dx < 2 / |VV^|Mx = 2 / IWlMx 

O JQiAsi) JAg, 

<8 I p'^\Vi/jfdx <16E^{u) <16A. (3.9) 
Ja„ 



The desired extension of ^ onto A is now given by ip{x) = ^(Q{x)). Using again 
the conformal invariance of the Dirichlet integral, we see by (3.9) that ^'^(e*^) = 
i Ia I^V^Pdx = ^ Jq \Vip\'^dx < 8A. Besides, since ip = ip on As/ and p = \u\ < 1, 
we get 



\u-e''^\\l2^A)= [ |M-e^^|Ma;+ /" {p-l^dx 

<A\A\As/\+ I {p^-lfdx<A{\A\A5,\+Ke'^). (3.10) 



As' 



Then, by choosing small 5' and Si — £i(5') > 0( £i(5') < s'l), in view of Lemma 10, 
bounds £;e(e^^) < 8A, E^{u) < A and (3.10) we have |abdeg(w) - abdeg(e^^)| < 1/2, 
when e < £i. But abdeg(e^^) = deg(e^^,aQ) = deg(e^^,>C) = deg(^,£), due to 
(3.2). Therefore if deg(^, C) ^ d then u e E^^'^, and vice versa. Thus (ii) is proved. 
□ 

4 Minimization among 5^- valued maps. Some up- 
per and lower bounds for problem (1.7) 

Consider the minimization problem 

Io{d, A') := mi{Eo{uy, u G H\A; S^), deg{u, duj') = deg{u, dQ') = d}, (4.1) 
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where Eq{u) = J^, \'Vu\'^dx, A' = Q' \ uj', and uj', Q' are smooth bounded simply 
connected domains in M^, such that u' C Q'. This problem is a particular case of 
the minimization problem considered in [14] (Chapter I). 

Proposition 13. [I4] There exists a unique (up to multiplication on constants with 
unit modulus) solution u of the minimization problem (4-1), o-nd u is a regular 
harmonic map in A' (i.e. —Ait = it|Vitp in A') satisfying u x ^ — Q on dA' . 

When A' = A, then any minimizer u of (4.1) belongs to J'dd- By (3.2) we also 
have abdeg(w) — d. This yields the following (optimal) bound for (1.7), in the case 
when p — q — d. 

Lemma 14. We have ms{d, d, d) < Io{d, A) for any e > 0. 

It is shown in [14] (Chapter I) that Io{d,A) can be expressed by 

Ioid,A) = ^J^\Vho\'dx, (4.2) 

where ho is the unique solution of the linear problem 

Aho = in A 

< Hq = 1 on dfl, ho — Const on dcu (4.3) 

/ ^da^2nd. 
KJdn diy 

ho and the solution V of (1.9) are related via ho = 1 + 27[d{V — l)/cap(A), where 
ceip{A) stands for iJ^-capacity of A (see, e.g., [25]). Thus Io{d,A) = 2{'7id)^/cap{A), 
and this clearly holds for any doubly connected domain A' in place of A. Therefore 
we have 

Lemma 15. Io{d,A') depends continuously on cap(A'). 

By using this simple result, we obtain the following lower bound for the GL 
energy of solutions u E J oi the equation (1.1). 

Lemma 16. There is £2 > such that for any solution u e Jim of GL equation 
(1.1) satisfying Es{u) < A, abdeg{u) E {d — 1/2, 0?+ 1/2) we have 

Ee{u) > Io{d, A)-^ + 7r{\d-l\ + \d- m|), (4.4) 

when e < 82, and £2 depends only on A. 
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Proof. Due to the maximum principle < 1 on A. As in Lemma 12 we consider 
the domain As> that is defined by (3.6) and depends on a positive parameter 6' < 1/2 
to be chosen later. Since |m| < 1 on ^4 we can apply Lemma 8 to get the bound 

\u\>l-e in As', (4.5) 

for £ < £2 (^2 — ^2{^': ^) > 0). Consider now the map 



1 \ u if ImI < 1 — £ 
u = { (4.6) 

(1 — s)]^ otherwise. 



1 -£ 



By (4.5) we have |-it| = 1 on As> and, according to Lemma 12, deg(-u, £) = d when 
e < min{£i,£2}- Consequently, the degree of u on both connected components of 
dA'g equals d, so that | J^^^ \Vu\'^dx > Io{d, As') (cf. (4.1)). Therefore, by using the 
obvious pointwise inequality \Vu\^ > 2|(?a;^'U x dx^ul and the integration by parts we 



get 



If If 



+ I / d^^ux dx^udx\> lQ{d,As') +Ti{\d-l\ + \d-m\), (4.7) 

k=l,2 ''^5' 

where A^^"^ and k — 1^2 are respectively the outer and the inner connected com- 
ponents of ^4 \ Ajt. On the other hand, it follows from (4.6) that \u\ < \u\ < 1. 
Therefore, 

E.{u) < -^^^^E,{u). (4.8) 

Bounds (4.7) and (4.8) yield (4.4) when 5' is such that /o(rf, Ay) > h^d^A) - 7r/4 
(cf. Lemma 15) and e is sufficiently small. □ 



5 From a vortexless minimizer to one with a 
single vortex 

The main Theorem 3 is proved by induction on the "number of vortices" in mini- 
mizers. More precisely, given integer d > 0, we show the existence of minimizers of 
(1.7) for p = g = d, then pass io p = d — 1, q = d and p = d, q = d — 1, e.t.c. The 
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key point of the proof is the induction step, when the degree changes by one on du 
or dVt. This change results in the rise of an additional vortex in a minimizer. For 
arbitrary p and g, satisfying the conditions of Theorem 3, this step is quite tech- 
nical. This is why we consider here a particular case of transition from p — q = d 
(no vortices) to p = d, q = d — 1 (one vortex). The transition from p — q — d to 
p — d — 1 and q — d is quite similar. We first establish 

Lemma 17. Given an integer d > 0. For sufficiently small e, e < with > 0, 
the infimum m^{d, d, d) in (1.7) is always attained, and me{d, d, d) < Io{d, A). 

Proof. Let u he a, weak if^-limit of a minimizing sequence {u^'^''). Since any 
minimizer v of problem (1.14) is an admissible testing map for problem (1.7), such 
a minimizing sequence exists and by using Price Lemma we obtain 

E,{u)+n{\l-d\ + \m-d\) <limmiE,{u^^'^) <- [ \Vv\'^dx = Io{d, A), (5.1) 

where / = deg{u,duj)), m = deg{u,dfl). Due to Proposition 2 we have abdeg(M) G 
{d — 1/2, d+ 1/2) when e < So, therefore the first variation of (1.2) at u vanishes, 
i.e. M is a solution of equation (1.1). Indeed, thanks to Lemma 10, we have for any 
w G Hq{A;M.'^) with sufficiently small iJ^-norm, u^''^ + w is an admissible testing 
map when k is large, hence Eg{u + w) — E^{u) = \imi^^Qo{Ee{u''^^ +w) — Ei;{u^^^) > 
(where lim^^oo denotes any partial limit), and we are done. Now, since -u is a 
solution of (1.1), we can apply Lemma 16. That is we substitute (4.4) in (5.1), this 
leads to 

\d-l\ + \d-m\ (5.2) 

when e < min{£o,£2}, i-e. I — m — d (since Z, m and d are integers). Thus the 
infimum in (1.7) for p = = d is always attained when e is sufficiently small. Lemma 
is proved. □ 

Next, we perform the transition from the minimization problem (1.7) for p — 
q — dto that ior p = d, q = d — 1 and show that me{d,d — 1, d) is always attained 
when e is sufficiently small. This is done by a comparison argument of ms{d,d—l, d) 
with the energy -^^(m) of a minimizer m of (1.7) for p = q ~ d. We first describe the 
properties of such a minimizer. 

In Section 5, it is shown that for small e any minimizer -u of (1.7) for p = g = is 
vortexless (see Remark 22), i.e. u = pe^^ with smooth p > and : /I — >• M \ 27r(iZ 
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(torus). It follows that we can write u = pe'"*^, the e*^ and V^^ being smooth maps 
globally on A. Then the boundary value problem (1. !)-(!. 3) rewritten in terms of p 
and 9 is 

'div(p2V^) = in A 



O] 

1 



— = on OA, 
ou 



-Ap + d"|V^|V+ ^P(P^ - 1) = in^ 



(5.4) 



p—1 on dA. 
We also have, in view of (2.2), 

V/i = -dp^W^e in A, (5.5) 

where h is the solution of (2.1). It follows that [h, 9) defines orthogonal local coor- 
dinates in a neighborhood of thus straightening out the boundary. Really, it is 
straightforward to verify that 

1 - h{du) = — \— I Vh ■ VVdx = 27rabdeg(ii)/cap(A), 

cap(yi) J 

while abdeg(u) > (d — 1/2) > 0. Then by applying the maximum principle to (2.4) 
we get, 1 > h{x) > h{dui) in A. This, in turn, implies, by Hopf's boundary lemma, 
that 1^ > on dQ; i.e. the map {h,9) : A ^ W X 'R\ 27rZ can be extended to a 
C^-diffeomorphism of a one sided neighborhood of dfl onto its image. Thus, there 
are some 5 > and a domain Gs C A such that 

xeGs^ih,9) e Us = {1-6,1) xR\ 27rZ 

is a one-to-one correspondence, which extends to a C^-diffeomorphsm of Gs onto 
[1 - 5, 1] X R \ 27rZ. 

The following Proposition is crucial for existence of minimizers of (1.7) for p — d, 
q — d — 1. In particular, combined with Lemma 14 it provides an independent of e 
bound for m^{d,d — 1, c?). 

Proposition 18. Let u = pc"'^, p > 0, be a minimizer of (1.7) for p = q = d. 
Assume that e is so small that Proposition 2 holds with A = Io{d, A). Then there is 
a testing map v e Jd{d-i) such that abdeg(M) e (ci — l/2,ci-|-l/2) and 

E,{v) - Ee{u) < TT. (5.6) 
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In Section 7, the generalized version of Proposition 18 is used to show existence 
of minimizers with several vortices. 

Proof of Proposition 18. We seek the testing map v in the form 

V = pwt, (5.7) 

with an unknown for the moment Wf. The following Lemma allows us to compare 
the energy E^{u) of u with that of v. 

Lemma 19. If w e H'^{G',R^), G' C A, is such that \w\ = 1 on G', then 

£(iv(p^)r + ^(Ip^p - 1)')^^ = ^/iv«r + ^(i«r - + 2Lf{w, go, 

where 

L'f^{w,G') = - [ p^\Vw\^dx-- [ |V^|W|'dx+-^ / p\\wf-lfdx (5.8) 
2 Jg' 2 Jq, Ae Jq, 

This result is a variant of the factorization argument due to [16], its proof is presented 
in the end of this section. 

Note that if G' — Gs we can rewrite the functional (5.8) by using local coordinates 
(h, e) as (cf. (5.5)) 

Lf\w,Gs)^'^ [ \d,w\'p'dhde 

- - (5.9) 

Instead of dealing with Li^\w, Gs) we will make use of the simplified functional with 
a quadratic penalty term, 

M^{w) = ^ I {d^\dhw\^ + \dew\^)dhde + ^ I (Alw-e^^-dVHdM^- (5-10) 
2" Jus 2a Jn^ 

This last functional admits the separation of variables. 

Now consider the map Wt that is given by Wt = c"''' in A\Gs, and continued 
to Gs as a minimizer of the functional M\{w), where A > 2rf^, with the following 
prescribed boundary data: 

Wt = e^''^:r,(e^^) on 9Q, (5.11) 
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wt = e"^ ondGs\dn, (5.12) 

where J^t{z) ■— Ct{z) (bar stands for the complex conjugate), Ct{z) — is the 

classical Mobius conformal map from the unit disk onto itself, i < 1 is a positive 
parameter. Both parameters A and t will be determined later. Since deg(^t, S^) — 
—1 and deg(e'^, dO) — 1, the standard properties of the topological degree implies 
that if V is as in (5.7), then 

deg(v, dVt) — d—1, deg(v, dto) — d. (5.13) 

The map Wt is well defined now, because the functional Mx{w) with Dirichlet condi- 
tion on the boundary has a unique minimizer for A > 2d'^. Moreover \wt\ < 2, since 
for if not then by taking Wt = ^ min{|u7(|, 2} in place of Wt the first term in (5.10) 
does not increase while the second one decreases, i.e. Mx{wt) < Mx{wt), it is a 
contradiction. 

Note that under the following choice of A, 

^^="^""{ 2.^infl|Vgp ''^l 
(|V^| > on the closure of Gs) we have 

p'{\wt\' - If < {\wt\ - lf{\wt\ + 1)' < \wt - e'^'lWwtl + If 

< 9\wt - e^'^Y < 2£'A|V^| V* - e'''T in Gs, 

thanks to the bounds \wt\ < 2 and p < 1. It follows that Li^\wt,Gs) < Mx{wt). 
We get now, by virtue of Lemma 19, that v — pwt satisfies 

E,{v)<E,{u) + Mx{wt). (5.14) 

We can obtain a representation for Mx{wt) in separated variables. Namely, 
expanding z'^J-'t{z) on as 

oo 

zdj^t(z) = (1 - t)z'' + t{t - 2) J](l - t)''z'^-''-\ 

k=0 

we have 

oo 

wt^il- tf^,{h)y + t{t - 2) ^(1 - tfh{h)e-^^'-<'+'^\ (5.15) 

A;=0 
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where fk{h) satisfy, according to (5.11, 5.12), 

/fe(l-5)=0, = (5.16) 

Substitute (5.15) into (5.10) to obtain 

J.2 J.2 °° 

M,{w,) = ^$-i(/-i) + ^ - 2)'(1 - tf'Mfk), (5.17) 

fe=0 

where 

$.(/.)= r (c^2|/^(/i)r + ((^-c^+l)2 + A-d2)|/,(/i)|2)d/i. (5.18) 
Minimizing (5.18) under the conditions (5.16) we get 

Qk+(h-l) gfc-Ci-l) 

^ 1 _ e(fc— fc+)<5 + i_e(fc+-*^-)<5' ^^"-^^^ 

where /c± = ±--^ (A; — d + 1)^ + A — d^. Therefore, we have 
a 

Mfk) = d{k-d+l){l + ^^ + 0{l/k^)), asA;^oo. (5.20) 
Finally, using (5.20) in (5.17), we obtain 

Mxiwt) < 7r((l - t)2 - If - + 27rt'(^ " J2 / + 

fe=0 fc=l 

= 7r(l - 2i - 2i'(A - d^) log(l - (1 - t)')) + (C + 7r)tl (5.21) 

Observe that the right hand side of (5.21) is strictly less than tt when i > is chosen 
sufficiently small. By (5.14), for such t the map v — pwt satisfies (5.6). 

It remains only to show that abdeg(^;) e(d — l/2,d+l/2). To this end note 
that by (5.15) and (5.19) Wt pointwise in as i — > 0. Therefore pwt — > 

pQide^_ weakly in H^{A), so that abdeg(pwt) — > abdeg(ii). On the other hand, 
by Proposition 2, we know d — 1/2 < abdeg(M) < d+ 1/2. Thus, after possibly 
passing to a smaller t, v = pwt satisfies the required property. □ 

Now we have that under the conditions of Proposition 18, there exists a minimi- 
zing sequence (u^''^) of admissible testing maps in problem (1.7) for p = d, q = d — 1 
such that lim^^oo E^{u^^^) < ms{d, d,d) + % and u^'^^ weakly i/^-converge to a map 
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u ^ J . Moreover, any minimizing sequence has a subsequence with the same pro- 
perties. Show that any weak hmit u is also an admissible map. Let / = deg(-u, Su;), 
m = dcg{u,dQ). By virtue of Lemma 16 (in the same way us in Lemma 17 one 
shows that u satisfies (1.1)) and Lemma 9 we have 

Io{d, A) - ^ + 7i{2\d - /| + |(i - 1 - m| + |(i - m|) 

< E^{u) + 7r(|ci - /| + Id - 1 - m|) < m,{d, d, d) + tt, (5.22) 

since abdeg(ii) = limfe_>oo abdeg(u*^'^)) e [d — 1/2, 1/2]. Due to Lemma 17 
m^{d, d, d) < Io{d, A) so that (5.22) implies that I — d and either m — d — 1 
or m — d. In the last case, u becomes an admissible map in problem (1.7) for 
p = q = d and therefore E^iu) > ms{d,d,d), which contradicts the last inequality 
in (5.22). Thus, u E Jd^d-i), abdeg(M) E [d — 1/2, d + 1/2], i.e. u is in the set of 
admissible testing maps of problem (1.7) for p = d, q ~ d — 1. 

Proof of Lemma 19. We have, by using (5.4), 

/ \V{pw)\Mx= I {p^\Vw\^ + Vp-V{p{\w\^ -l)) + \Vp\^)dx 

JGs Jgs 

= / (p^iwwi' + d'p'iwei' + ^p'ip'-mx 

- 1 (dviveivi'+ V(p'-i)kr-ivpr)dx. 

JGs ^ 

Then simple algebraic manipulations give the required result. □ 



6 Asymptotic behavior of local minimizers 

In the previous section, we established the existence of minimizers of (1.2) in jTj^-' 
and demonstrated the first induction step of the proof of Theorem 3 that consists 
in transition from p — q — d io p — d, q — d — 1 m. {1.7). (In fact, modulo the 
assumption that minimizers in J7j^^ are vortexless, we actually proved the existence 
of minimizers in J^^fLi^d ^t^d-iy) ■'■^ order to show the induction step for any 
integer p < d and q < d,we need to establish some properties of minimizers of (1.7), 
and we are especially interested in their behavior near the boundary. At this point, 
we assume we are given a family {u^} of minimizers for (1.7) and 

E,{u,) < A := Io{d, A)+7r{\d-p\ + \d-q\). (6.1) 
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Also, wc suppose also e < Eq, where Sq = eo{A) > is as in Proposition 2. It follows 
that maps are local minimizers of Es{u) in J' and therefore they satisfy (1.1), 
(1.3). 

We will use the notations: -Bed/) = ^ • \^ ~ vl < Pe{x) = \'<J'e{x)\, 
he{x) is the unique solution of (2.1) (associated to m^), and C is the contour as in 
Lemma 12. The contour C separates the two open subdomains in A, where 
is the domain enclosed by dVl and C and Q~ — A \ [Q^ U C). We also set 
Qf^{xeQ^;pl{x)<l-e^'-'}. 

6.1 Proof of Theorem 4 

Since |V/ie| < IVm^I (by Lemma 6), the family {he} is bounded in H^[A), and 
therefore there is a sequence — > such that 

he^^ — >■ h weakly in H^[A), as A; — > oo. (6.2) 

In order to identify /i, we make use of Lemma 8 to obtain that, up to a subsequence, 
maps u^^ converge to a yS^- valued map u in C^^J^A). Since d^^u x d^^u — a.e. in 
A, we have ^h^^ — 2dxj^u^^. x d^^u^^ ^ in C^^J^A), thus /i is a harmonic function. 
Moreover, h — 1 on dO, and h — Const on du. On the other hand, 

abdeg(KeJ = 7^ / V/ie, • VVdx ^ — / V/i • VV^da; = — / ^ds. 

According to the property c) of abdeg( ■ ) (see Proposition 11 in Section 3), 
abdeg(M£) — >• d, as £ ^ 0, therefore h = Hq (where ho is the unique solution of 
(4.3)) and the convergence in (6.2) holds for the whole family {he}- Thus, applying 
again Lemma 8, we obtain 

he — >• ho in Cioj,(A) and weakly in H^{A), as £ — >• 0. (6.3) 

By (6.1) and Lemma 8, maps converge, up to a subsequence, to u & H^{A; S^) 
in Cl^^{A) and weakly in H^{A). Moreover, abdeg(M) ~ d and in view of (6.3) 
|Vm| = |V/io| a.e. in A. It follows that is a solution of the minimization problem 
(1.14). 

In order to demonstrate the energy convergence stated in Theorem 4, we argue 
as follows: by using two pointwise equalities | Vwep = 2dx{ae x 8x2^^ + 4|(92M£p and 
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iVwep = —2dx^Ue X Ox^Ue + 4|c?2M£p and the pointwise inequality \Vue\ > jV/iel, we 
have 

^ / \Vue\^dx > - / dx^Us X dx^Uedx + 2 / \dzUs\^dx 

^ J A JQt JQt 

+ / dx^u^ X dx^u^dx ^ 2 I \d^u^\'^dx + I- [ \S/h^\'^dx. (6.4) 

Jqi Jq7 2 Ja\{q+vjQ7) 

Let us estimate the right hand side of (6.4) from below. Introducing ae{x) — 
max{p2(x), 1 - £1/2}^ we have (by (2.3) (2.4)) 

rv ( 1 v/.^ 2 [O in A\{QtUQ-)- 

div(— -V/iJ = YZTTPI \ (^•^) 
'^^v^/ ^ I dx^Ue X 9j;2it£ otherwise. 



Integrating (6.5) over we get for sufficiently small £, 



^-ri^ / 5a;i«£ X dx^Uedx = f ^ds - f 

1 - Jot Jan ou Jc 



dhe ds 
c dv pI{x) 

Ue X ^— ds - -. — r X —- — -ds = 27r{q - d), 
an Jc Pel Pel 



where we have used Lemma 8 and Lemma 12. Thus, we have 

/ dx^u^ X dx2U^dx = (1 - £^/^)7r(g - d). (6.6) 
JQt 

Similarly, integrating (6.5) over Q~ we obtain 

/ dx.ue X dx.uedx = (1 - e'/^)7r{d - p). (6.7) 
In order to estimate the last term in the right hand side of (6.4), we write it as 
/ \Vh^\^dx = / \Vh^ - V/iol^dx 

JA\{Q+liQ7) JA\{Q+uQe) 

+ / {2Vhe - Vho) ■ Vhodx - / {2Vhe - Vho) ■ Vhodx, 

J A JQtLlQe 

and note that by virtue of (6.1) the measure of Qf U Q~ vanishes as £ ^ 0, so that 
/ \Vhe\^dx= I \Vh,-Vho\'^dx+ \Vho\'^dx + o{l). (6.8) 

JA\{Q+LlQe) JA\{Q+VjQe) Ja 

Thus (6.6-6.8, 6.4, 6.1) imply E^{ue) Eq{u) + 7r(|d - p| + |rf - g|). □ 
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As a byproduct of the above proof by, (6.6)-(6.8), (6.4) and (6.1) we get 

/ {\u,\^-lfdx^o{e'), (6.9) 

J A 

I |V/i£- V/io|^dx = o(l), (6.10) 

JA\{Q+yjQ7) 

I \dzUe\^<ix ^ o{l) I \d-zUe\^dx^o{l). (6.11) 

Jot JQ7 



6.2 Properties of minimizers of (1.7) for small e 

First, by using (6.9) and the following methods of [14] we get that converges to 
1 uniformly on compacts in A. Moreover, we have 

Lemma 20. For any fx > we have 

sup{dist(y, dA);ye A, pl{y) < 1 - /x} = o{e). (6.12) 



Proof. Assume by contradiction, that for a sequence £fe — > and 7 > we have 
pl^{yk) < 1 - A* and dist(yfe,9A) > 7^^. Due to (2.6), iVjiieJ^I < a/sk in BxekiVk), 
where < A < 7 and a{— a (A)) is independent of Sk- It follows that \ue^{x)\'^ < 
1 — 11 + 5a when x e Bs^^{yk) and 5 < \. Then Bsei^{yk) C A and 

\ [ iKJ^ - ifdx > n{p - SafS^ > 0, 

as soon as < 5 < min{A, ///(2a)}. This contradicts (6.9). □ 
Important properties of and h^, in a vicinity of the boudary dA, are established 

in 

Lemma 21. For any < p, < 1 and k < 1 there are ei{p),e2{p, k) > such that if 
pHv) < 1 ~ then 

(a) fore < ii{p) we have hs{y) > l + p/i if dist{y,dfl) < e andh^{y) < r^ldcu) — 
p/A if dist{y , duj) < e; 



(b) for e < £2 (a*, i^) we have 

1 

2 



\Vu,\'^dx > KTT. (6.13) 
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Proof, (by contradiction) Let us assume that either (a) or (b) is violated for a 
sequence St ^ and some y = Uk such that hf.^{yk) < 1 — /U. According to Lemma 
20, yk —>■ OA. For the definiteness we suppose that yk —>■ dfl, then (by Lemma 20) 

dist{yk,dQ) ^o{ek). (6.14) 

Let be continued in u in such a way that < CH-UgH and l-u^l < 

1 in Q, where C is independed of e. We also assume that = h^{duj) on uo. 
Following [12] we rescale m^j, and /i^^ by a conformal map that 'moves' yk away 
from the boundary. Fix a conformal mapping rj from VL onto the unit disk i?i(0). 
We introduce the conformal map C,k{z) = (-2 — fliUk))/ {'n{yk)z — 1) from -Bi(O) onto 
itself and set Uk{z) = Ue^,{i]'^{Ck{z))), Hk{z) = he^{r]~^{Ck{z))). It is easy to see 
that ||t/jfc||i/i(a) < C and ||i?jt||iji(n) < C with some C independent of k. Therefore, 
without loss of generality, we can assume that Uk and if ^-weakly converge to 
limits U and H, respectively, as /c — > oo. 

Arguing, as in [12] (Section 4), we can show that Uk ^ U in Cl^^{Bi{0)) and that 
AU = OinSi(O). Therefore, |C/(0)|2 = lim^^oo |C/ik(0)|' = ^^k-.oo MVk)]'' < 1-/^- 
We also have \U\ — 1 a.e. on S^. Show now that dzU = in -Bi(O). Indeed, by 
the maximum principle \U\ < 1 in Bi(0), hence maxBt{o) \Ukiz)\'^ < 1 — e]/"^ for 
any fixed < i < 1 and sufficiently large k. It follows that for such k we have 
V~^{Ck{Bt{0))) C Qf. Then in view of (6.11) we get, by using the conformality of 
the maps r}~^ and Cfc, 



k\ ax = / l^zUckl / l^zUck 

B*(0) Jv-HUBtm) J At 

This implies that d^U = in fii(O). 

In order to show (b) we use the pointwise equalities ||V^7p = —dx^U x dx^U + 
and dzU = to obtain 

]- f \VU\Mx = - /" dx,U X dx^Udx = -7rdeg(C/, S^). 

As U ^ Const, we therefore have ^ /g^^Q) |V[/|^da; > tt. It follows that there is 
< i < 1 such that 

- / \VU\'^dx> KTT. (6.15) 

The image Cfe(i?f(0)) of the disk Bt{0) is the disk Bt^{^k) with the radius tk = 
^SiSjp center at ^k = • According to (6.14) tk = o{sk) for 
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k oo, hence ri^^{Bt^{^k)) C Bi;^{yk) when k is sufficiently large. Then, by using 
the conformal invariance and lower semicontinuity of the Dirichlet integral, and 
bound (6.15), we get 

/ iVwejMx > / |Vw£^|Mx= / |VC/fe|Mx > 2K7r as A; ^ oo. 

In order to show that hs^{yk) — Hk{0) > 1 + /x/4 when A; — > oo we note that the 
system (2.1) is conformally invariant, i.e. 

V^Hk = {UkXd,,Uk,UkXdM in Ck\v{A)). 

Then, bearing in mind the convergence properties of Uk, we obtain that if^ — > if in 
CUBim and 

V^H={Uxd,,U,Uxd,,U) = ~V^{\U\^) in B^{0), 

where we have used the fact that d^U = in -Bi(O). Since H — \U\ — 1 on dBi{0) 
we have H — | — in Si(0), therefore 

hm hM = hm H,{0) = I - l\Um' > 1 + ^. 

fe— ♦oo K— »oo L L Z 

Lemma is proved. □ 

Remark 22. Lemma 21 implies that in the case when p = q = d minimizers of 
(1.7) are vortexless for sufficiently small e. Really, by Theorem 4 they if^-strongly 
converge, up to a subsequence, as £ —> to a minimizing harmonic map u e ■ 
On the other hand (6.13), exhibits the energy concentration property near zeros of 
minimizers, which is incompatible with the strong if^-convergence. 

The following result, describing the structure of the function he for small e plays 
a crucial role in the proof of the main technical result (Lemma 25) in Section 7. 

Lemma 23. We have, for small e, e < €4 (where £4 > Oj, 

(i) pI{x) > 1/2 when he{du) - 1/8 < h^ix) < 9/8 and K{du) < miucKix) < 
m.ax.che{x) < he{dfl), while if p1{x) < 1/2 then either 

d\si{x,dQ) < dist(£,aO) and h^ix) > 9/8 

or 

dist{x,du) < dist(>C, 9a;) and he{x) < hs{duj) — 1/8; 



29 



(n) there are a;* G dVt, a;** G duj such that ^{xD > and f^(a;**) > 0. 

Proof, (i) follows from Lemma 20, Lemma 21 and the convergence properties of 
as £ — > established in the course of the proving Theorem 4. To prove (ii) we argue 
as follows. Let £2(a*, be the best (biggest) constant in Lemma 21. Then £2(a*, 
is increasing in // and decreasing in k. For A; = 1, 2, . . . set 

= when min{£2(fe^, fc^),l/(A; + l)} < £ < min{£2(|,¥)'lA}- 

Then /i^ — ^ as £ — >• and (6.13) is satisfied with n = 1 — (1^ when pliy) < 1 — fie', 
the same being true when /ig is replaced by — max{/i£, e^^^}. We pick a point xi^^ 
in A such that p1{x^P) < 1 — jie'-, then we pick a point xf^ in A \ B2£{x'i^) such that 
p1{x'i^) < He , e.t.c. unless for some we have p^ix) > l — pe on A\\J^^-^^B2e{xf^). 
By the construction of since disks Be{x^e^) are disjoint, 

^ / \Vu'\'^dx > ^ V / IVli'l^dx > Ke{l - Pe)T^. 

Therefore by (6.1) we have a uniform bound < C. Arguing as in [14] (Chapter 
IV, Theorem IV. 1) we can increase the radii of disks to e\ > 2e (with A independent 
of e) and take a subset of {1, ... , K^} such that 

\JkehBex{.xl) D uf^iB2e(x^) and dist( x^ ,x^) > 4£A for different k,k' G I^- 

We also have 



pi{x) >l-peOnA\ {JkaeBsxixD 

Assume and extend to the whole and set — he{du), Hq — ho{du;) 
in u, and ^ ho ^ I in \ Q. Since pl{x) > I - Pe > ^ - e^^^ in ^ = 
B2xM''^) \ Bxe{xf^). by (6.10) we have 

/ \Vhe\^dx<2 (|V(/ie -/io)r + iV/ioDdx ^ as £ ^ 0. 

Then, writing the integral over D^"^ in the polar coordinates with the center at x^e"^ 
and using Fubini's theorem, we can find Xf\ Xe < < 2Xs, such that 



\Vhe\^da = o{l/s) 
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Therefore, by the Cauchy-Shwartz inequahty, 

/ |V/i,|d5<(7rAf))V2| /■ jVh^lMsY' = o{l). (6.16) 

Now, integrate(2.4) over Q^\Uk^i^B^(k){xi''^) to get, according to (6.16) and Lemma 
12, 

/ ^ds = 27rd + > / = 27rd - o(l) when e ^ 0, 

Jr. t^J\x-xi'^\=xi'' du 

where I'^ denotes the subset of indexes k e such that B^(k){xi''^) fl Q'^ ^ 0, and 
= \ \Jk(ivB^(k){xf^). Therefore, there is x% e such that ^{x%) > 0. 
Similarly, we can show that on — doi)\ Ujtg/^\//-B^{fc) (xi*^^) there is x** such that 

^{xT) > 0. ' □ 



7 Inductive proof of Theorem 3 

Fix A > and an integer d > such that A > Io{d,A), and let aeo be the greatest 
integer such that 

Io{d,A) + Traeo < A. 

Clearly, aeo > 0. In this Section we show that for small e the infimum in problem 
(1.7) is always attained, provided integers p, q satisfy 

36(p, ?) < aeo and p < d, q < d, (7.1) 

where ae(p, q) — \d — q\ + \d — p\. 

Proposition 24. Given an integer K < ^o- Let p < q, q < d he integers such that 
^{p,q) < K. Then, for sufficiently small e the infimum in problem (l-l) is always 
attained and 

me{p,q,d) < m£{l,m,d) + 'K{{l—p) + {m — q)) when p<l<d,q<m<d. (7.2) 
Moreover, the inequality in (7.2) is strict unless I — p and m — q. 

Proof. The proof is by induction on K. The basis of induction (K=0) is established 
in Section 4 (cf. Lemma 17). The demonstration of the induction step relies on the 
following Lemma, whose proof is in the end of this section. 
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Lemma 25. Assume integers p and q satisfy (7.1), and for e < e^, £5 > 0, there 
exists a minimizer of problem (1.7) whose GL energy E^i^u^) satisfies the bound 
(6.1). Then for any e < m\n{ei,e^} (where 64 is as in Lemma 23) there exists 
Ve G ^(q-i) such that abdeg{vs) e {d — l/2,d + 1/2) and 

Ee{ve) <m,{p,q,d) + n. (7.3) 

Similarly, in J^(p-i)q there exists a testing map (still denoted v^) satisfying (7.3), and 
such that abdeg(f£) e (d - 1/2, d + 1/2). 

In view of Lemma 25, to prove the claim of Proposition 24 loi K + 1 in place of 
K it suffices to show that mg(p,q — l,d) is always attained for sufficiently small e 
whenever p < d. q < d, and si{p, q) = K (the attainability of m^{p— 1, g, d) is proved 
similarly). Let u be a weak i^f^-limit of a minimizing sequence {u^'^^). According 
to Lemma 25, such a minimizing sequence exists, moreover, this lemma, with the 
induction hypothesis, imply 

limsup E,{u^^^) < rueil', m', d) + 7r(/' - p) + 7r{m' - {q - 1)), (7.4) 
fe— +00 

where p<l'<d, q — l<m'<d and ee{l', m') < K. We know that abdeg(ii) = 
limfe_>oo abdeg(ii^'^^) e (d — 1/2, d + 1/2) when £ < £0 (where £0 = £o(A) is as in 
Proposition 2), hence is a solution of the GL equation (1.1) (see arguments in 
Lemma 17). Therefore, if we write \\m.'m.lk^ooE^{u''^^) < \\m.suY)y._^^Es{u^^'>) and 
apply successively Lemma 9 and Lemma 16 to the left hand side, we get by using 
(7.4) with I' — m' — d that for e sufficiently small 

TT 

Io{d,A) - - + nse{l,m) + n{\l - p| + |m - (g - 1)|) < m^{d, d, d) + 7rae(p, q — 1), 

where / = deg(M, 5a;), m = deg(M, dQ). Thanks to Lemma 14 me{d, d, d) < Io{d, A), 
thus 

\l-d\ + \l-p\< |p-d| + l/2and \m - d\ + \m - {q - 1)\ < |(g - 1) - d| + 1/2. 

Since / and m are integers, it follows that p<l<d, q— l<m<d. Now, assuming 
/ 7^ p or m 7^ g — 1, we use Lemma 9 and (7.4) with I' = I, m' — m to obtain the 
following 

Ee{u) + 7r{l -p) < liminf ^^(m'^''^) < me{l,m,d) + n{l - p) + n{m - (g - 1)). 

k—*oo 
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On the other hand, u G jTjm and abdeg(-u) G [rf — 1/2,6? + 1/2]. Hence Es{u) > 
mir{l,m,d), which is a contradiction. Therefore / = p and m = q — 1, i.e. u is 
an admissible testing map in problem (1.7) and thus the infimum me{p, q — l,d) is 
always attained. □ 
Proof of Lemma 25. For simplicity we drop subscript e. The underlying idea is 
to modify the minimizer u of (1.7) in a neighborhood of dfl as in Proposition 18 
(see Section 5). In general u is with zeros now, thus the arguments need to be more 
sophisticated. Loosely speaking, we construct a testing map v with an additional 
"vortex" located "near" x*, where x* is a point on dQ such that ^{x*) > (cf. 
Lemma 23). 

Step 1: Domain decomposition. Let 1 — 5, where 5 > 0, be a regular value of 
h (thanks to Stard's lemma this holds for almost all S). Consider the subdomain of 
A where h > 1 — 5. There is a (unique) connected component Ds of this subdomain, 
such that dDg D dfl. Since h{dQ) = 1 > h{duj), when 5 is sufficiently small the 
boundary of Dg contains a connected component ^ dQ. enclosing uj. According 
to Lemma 23, we can choose 5 small enough 5 < 5q {5q > 0) so that the domain 
enclosed by F^ and dO, hes away from the contour C, i.e. F^ also encloses C, moreover 
if p{x){= \u{x)\) vanishes at a point x of this domain, then h{x) > 1. Therefore, the 
minimum of h over the closure of the forementioned domain cannot be attained at 
any interior point, otherwise h satisfies d\y{-^h) = in a neighborhood of this point, 
which is impossible. In other words, h > 1 — 5 in the domain enclosed by F^ and dQ, 
i.e. this domain coincides with Ds- Thus, the boundary of Ds consists of exactly 
two connected components dQ and F^. Also, possibly choosing smaller 5q, we have 
that the set P = {a; G Ds; h{x) > 1} is independent of 5 (recall that 5 < 5q and 
1 — 5 is a regular value of h). Indeed, consider the set Ss = {x ^ A; h{x) > a} r\Ds 
(0 < 5 < ^o); where « is a regular value of h and 1 < a < 9/8. Ss consists of a 
finite number n{5) of connected components. Since Ds D Ds' ii 5 > 5', the function 
n{5) is nondecreasing, hence n{5) = \ims'^on{5') when < 5 < 5q (for some 5q > 0) 
and Ss = Ss' if 5,5' G (0,5o)- It follows that 1 — 5 <h<a in Ds\ Ds' when 
< 5' < 5 < 5o. For such 5 and 5' the function h satisfies div(^/i) = in \ Ds' 
(by Lemma 23) while h < 1 on the boundary of Ds \ Ds' , hence h < 1 in \ Ds' . We 
see now that {x G Ds; h{x) > 1} — P ■— r\s'<SoDs', when < 5 < Sq, as required. 
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Thus we have, for 6 < 60 

,, 1 f dhds 1 f u d u , , , ^ , 

d' :=— ^ = 7r ri^ TT^^* = ^^S^' > 

zvr J-p^ OP zvr J^^ \u\ or \u\ 

(1 — (5 is a regular value of h and h> 1 — 5 in D^) and the integer d! is independent 
of 5. Therefore u admits the representation u = pe^'^'^ in Gg = Ds \ P, where 
6 : Gs ^ M. \ 2ttZ is a smooth function and p > in Gs- 

Without loss of generality we can assume that u{x*) = 1. Since Vh = —d'p'^'V^O 
in Gs and §^(a;*) > 0, the map x {h,9) from a neighborhood G'^ of x* onto its 
image h{G'^) is a C^-difFemorphism. Choosing 5 small enough, we can assume that 

is defined by 

xeG'^ X eGs, 1-5 < h{x) < 1, e{x) E (-5,(5)(mod27rZ). 

Now we have A = G'g U G'^ U {A\ Gs) (see Fig. 1), where G'^ = Gs\ G'g. 




Figure 1: Domain decomposition. 

Step 2: Construction of the testing map. We seek the tasting map v in 
the form 

{u in A \ Gs, 
(7.5) 
pwt m Gs, 

with (unknown for the moment) Wt = Wt{h{x), 6{x)). Impose the following boundary 
conditions on dGs, 

wt = e''''' on r^, (7.7) 
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where < (/? < 1 is a smooth 27r-periodic cut-off function such that ^p{9) = 1, when 
9 e {-6/2,6/2) (mod 27rZ) and ^{9) = if ^ ^ (-5, 5) + 27rZ. It is easy to see 
that if Wt (considered as a function of h, 9) satisfies (7.6) and (7.7), when h = 1 and 
h = 1 — 6, respectively, and is a smooth 27r— periodic in 9 map defined in the strip 
1-6 < h < 1 then (7.5) defines for < i < 1 a map v e H'^{A; R^) such that 1^1 = 1 
on dA and 

deg(v, do,) — q — 1, deg(v, dcv) — p. (7.8) 
Expand the right hand side of (7.6) into the series 

e'^'-{l-M9)) 



^ = (1 - tb^i{t))e''^'^ + 1 J2 h{t)e-'^^-'^'+^'^\ (7.9) 
and set 



wt{h,9) = {l-tb_,{t)f_,{h))e^''' + t J2 h{t)h{h)e-^^'-''^'^', (7.10) 



where functions are defined by (5.19) with k± — \J {k — d' -\- ly -\- X — d'^. 
The positive parameters t < 1 and A > 2d'^ are to be specified later on. In what 
concerns the coefficients bk{t), we have 

\bk{t) - Ck{t)\ < C(l + |A;|)-^ Vn > 0, (7.11) 

where C = C{n) is independent of t, and Ck = {t — 2)(1 — tY for A; > 0, c_i = 1, 
Cfe = for A-; < — 1. The estimate (7.11) is obtained in a standard way, by comparing 
the Fourier coefficients in (7.9) with those of c'^'\e-i<^ - (1 - t))/(c"^^(l - t) - 1). 

Step 3: Verification of (7.3). Thanks to Lemma 19 we have (since \wt\ — 1 
on dGs due to (7.6, 7.7)) 



E,{v)=E,{u)+Lf\wt,Gs) 



where the functional L^f \w) is defined as in (5.8). Let us show, that for sufficiently 
small t 

Lf\wt,Gs) <7r. (7.12) 
To this end, like in the proof of Proposition 18, consider the quadratic functional 



liwt) = J / {d"\drwt\' + \doWt\' + X\wt- e^T - d''\wMV9\'dx. (7.13) 



'Gs 

35 



(One can actually show that Wt minimizes functional (7.13) under the boundary 
conditions (7.6), (7.7).) We have, since V/i = —d'p^V^O in Gs and p < 1 in A, 

f p^\Vwt\Mx< I {d'^\drWt\^ + \deWt\^)p'\Ve\^dx. (7.14) 

Jgs 

Moreover, if we put 

A = max \ — — . ^ , 2d'^ 

I 2£^mmg^ \ \/0\^ 

under the additional assumption that \wt\ < 2 in G^, the following pointwise inequa- 
lity 2e'^X\wt — e'^p|V^p > p^{\wt\'^ — 1)^ holds in G'g (see the proof of Proposition 
18). Thus, 

To demonstrate (7.12) we first note that 

J.2 °° 

= ^ E Mt)\'^'kifk), (7.16) 

fc=— oo 

where functionals are defined as functionals $fc in (5.18) with d' in place of d. 
The representation (7.16) is analogous to that in (5.17) while its justification differs 
because of the fact that V/i vanishes at least at some points of the boundary of Gg 
(and possibly somewhere in Gg). We rely on the following Lemma, which implies 
directly (7.16). 

Lemma 26. Let f,gEC^{[l — 6,l]; C) then for all integers n, m 

p 0, if n m 

J /(/i)e'"^^(/i)e^™V'|V^|Mx= Ktt 



'Gs \~d/^ I f(^^9{s)ds, if n = m. 
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Proof. By virtue of the pointwise equalities Vh ■ Vd = and div(p^^) = in Gg, 
for any regular values a, [3 oih such that \ — 5<a<l3<l and any integer 
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/e^"^^e^"^V'|V^|Mx = — [ Ve ■ vS^^^fgp'dx 

a<h<p n — m Ja<h</3 

^ ^ div(p'Ve)/^e^("-"^)^dx 



n — m 



a<h<(3 

+ I Ve -Vhif'g + fg')Q'^"'^^ p'^dx ^Q, (7.17) 

n — m Ja<h<p 
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where all integrals are understood over subsets of Gs- If n = m we set F{h) — 
£f{s)g{s)ds then, since \Vh\ = d'p'^\V9\ 

[ mmp'Mi^dx = ' [ viFih)).vh^='FiP)[ §^ 

Ja<h<p d' Ja<h<(3 P d' Jh=f3 P 

- 1^ / divi-Vh)Fih)dx = - / /(.)^(.)d., (7.18) 

here we have also used the fact that div(^V/i) = in G5. The statement of the 
Lemma is then obtained by passing to the limits a — > 1 — 5 and /3 — > 1 in (7.17, 
7.18). □ 
By using (7.11) in (7.16) we compute 

2 00 00 

^aK) = ^ E \ck{t)\'mk)+o{t')=7r{{i-tr-irj2k{^-tr 

fc=— 00 k=0 

00 +\2k 

+ 27Tt''{X - d'^) E ~, + 0{f) ^n-2nt + o{t). (7.19) 

k=l 

Explicit, but tedious computations (left to the reader), show also that 

\wt - e'*^'^! < Ct when ^ (-a, a) + 27rZ, (7.20) 

for any < a < 27r, where C is independent of t. Prom (7.20), we see the second 
term in (7.15) is of order 0{t^) as t ^ 0. Combined with (7.19) this proves (7.12). 

Final step. The bound (7.15) is shown assuming that l^tl < 2 in G'^. Note, 
that this can always be achieved by replacing wt by Wt '■— 'J^'t min{l, 2/1^^1}, and 
this change increases neither the first term in (7.15) nor the second one. Thus in 
order to complete the proof of the lemma, we need to show only that the map v 
defined by (7.5) satisfies d — 1/2 < abdeg(v) < d+1/2 when t is chosen sufficiently 
small. Indeed, due to (7.20) Wt weakly iJ^-converges to e"^'^. Therefore the norm 
\\u — tends to when i — > 0. Then, according to Lemma 10, for small t 

abdeg(v) is close to abdeg(M), while d — 1/2 < abdeg(M) < d+1/2 and we are done. 
□ 
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Appendix A 

Here is a simple example of Mobius (Blashke) test map that illustrates an important 
property of nearboundary vortices. Namely, each such vortex appears in pair with 
a "ghost" antivortex (vortex with opposite sign) coming from outside the domain. 
For simplicity consider the minimization problem of (1.2) with A — Bi, boundary 
condition = 1 on dBi and imposed degree one on the boundary dBi. It is shown 
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in [11] that the infimum in this problem is not attained and the behavior of the 
minimizing sequence is described by 

\Ce\CeZ-l \Cs\ Z - iCe + 0{l)y 

where Q G -^i, ICel ^ 1 as £ ^ 0. We assume that dist(C£, dBi) — o{e). Clearly this 
map has single zero of degree one at z = Ce- Introduce (* — 1/Ce, and write this 
map as 



Ve{z) 



1 


\z-C\ 


Z-Ce 




\z-Q\ 


\Z-Ce\ 



i^e{z) 



\Z-Ce 



CI 

Q\ 



(A.l) 



where u!e{z) is a real-valued function. Then, the last factor in (A.l) corresponds 
to a "ghost" vortex with the center at Q ^ Bi. The complex conjugation means 
this vortex has degree —1 which is why it is called antivortex. Finally, uJs{z) — 
1 + o(l) outside the disk of radius e with center z = Ce- These observations show 
superposition of vortex and antivortex has almost no energy away from the core, 
whereas, the contribution in the energy from outside the core of an inner vortex is 
known to grow as log(l/£). 



Appendix B 

Let us show that, for any integers p, q, d, the set of admissible testing maps in 
problem (1.7) is not empty and 



me{p, q, d) < loi^d, A) + n\q — d\ + n\p — d\. 



To this end chose two sequences {xf^^), (x^^) C A such that x^^^ dfl, x 
ier solutions pi'\ p^^ oi 

Apf^ = 2n\d - q\5{x - x[^'^) in A 



(fe) 

2 



(B.l) 
do;. 



(k) (k) 

Consider solutions Pi , P2 of the boundary value problems 



pf ^ = on dn, 



Pi 



(k) 



Const on du 



and < 



dp'; 



(k) 



-da = 2n\d — q\, 



Ap[^'^ = 2n\d- p\6{x - xi^^) in A 
p^^ = on duj, 

Const on dQ 



(fe) 

P2 



du 



da = 2n\d — p\. 



an du 

One can show that, for any neighborhoods U and V of dQ and dco in A, 

in C\A\U) and p^''^ ^ in C^(^ \ l^). 



(ifc) 
Pi 



(B.2) 
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There exist harmonic conjugates 

:A\{xP}^R\27r\d-q\Z, 
V'f :A\{4'^}^M\27r|d-p|Z, 

(i.e. Vt/jj'^ — V^pf \ j = 1,2) so that e''^ ^+'>''^\ ^ are harmonic functions in A. By 
using the pointwise equahty ||Viip = d^^ux dx2U-\-\\dzu\'^ and integrating by parts 
we get 

]- [ e^^i'^lVpf l^ + lVV^f^ndx^^ / iVe^i'^+^'^i'^dx^TrM-gl. (B.3) 
2 ^ Ja 

Similarly we have, 

lJ^e'^'^\\Vp^^^\' + \v4'^\')dx^n\d-p\. (B.4) 
Then we introduce 

where -u is a minimizer of problem (1.14) and V is the solution of (1.9). From (B.2- 
B.4) we derive that, for k sufficiently large, ■u*^'^^ is an admissible testing map for 
problem (1.7), and that E^{uP^^) — > /o((i. A) + 7r|g — d\+ tt\p — d\ when k ^ oo. 
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